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We prove that the stochastic and standard field-theoretical approaches produce exactly the same 
results for the amount of light massive scalar field fluctuations generated during inflation in the 
leading order of the slow-roll approximation. This is true both in the case for which this field is a 
test one and inflation is driven by another field, and the case for which the field plays the role of 
inflaton itself. In the latter case, in order to calculate the mean square of the gauge-invariant metric 
fluctuations, the logarithm of the scale factor a has be used as the time variable in the Fokker-Planck 
equation in the stochastic approach. The implications of particle production during inflation for the 
second stage of inflation and for the moduli problem are also discussed. The case of a massless self- 
interacting test scalar field in de Sitter background with a zero initial renormalized mean square is 
also considered in order to show how the stochastic approach can easily produce results corresponding 
to diagrams with an arbitrary number of scalar field loops in the field-theoretical approach (explicit 
results up to 4 loops included are presented). 

PACS numbers: 04.62,+v, 98.80.Cq 



I. INTRODUCTION 

It has been known for a long time that light minimally 
coupled scalar fields typically have anomalously high vac- 
uum expectation values (VEV) for even powers of fields in 
the de Sitter background. In particular, the mean square 
of a free, massive and minimally coupled to gravity scalar 
field in the equilibrium, de Sitter invariant quantum state 
(the Bunch-Davies vacuum) is [l[ 



3H 4 
8n 2 m' 



> H 2 



(1) 



if m 2 <C H 2 . Here H = a /a is the Hubble parameter 
and a(t) is the scale factor of a Friedmann-Robertson- 
Walker (FRW) cosmological model. For the particular 
case m = and (</> 2 ) = (or small) at the beginning of 
the de Sitter stage (set at t = here), we have 
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which grows without bound @, H, 0]- This anomalous 
growth is an essentially infrared effect, it occurs due to 
field Fourier modes with wavelengths far exceeding the de 
Sitter event horizon. Thus, it is not a consequence of the 
Gibbons-Hawking effective temperature T = H/2n [|[ 
experienced by a point observer inside her/his de Sitter 
horizon. In particular, in contrast with temperature- like 
effects which should be universal for fields of all spins, 
this effect occurs for minimally coupled light scalar fields 
and gravitons only (but not for, e.g., photons or light 
fermions). 



Because of continuous creation of infrared modes and 
growth of their occupation number, the quantum scalar 
field can be split into a long wave (coarse grained) com- 
ponent and a short-wave (perturbative) one. Then it can 
be proved that the former component effectively becomes 
quasi-classical, though random (i.e. all non-commutative 
parts of it may be neglected), and it experiences a ran- 
dom walk described by the stochastic inflation approach. 
In some specific case but beyond one-loop approximation, 
this approach was used already in ; see 6| for the rig- 
orous derivation in the generic case when this scalar field 
is a slow-rolling inflaton itself and for a number of ana- 
lytic non-perturbative results.* The case of a scalar field 
with the quartic self-interaction in the exact de Sitter 
space-time was first studied using this approach in 0], 
and a number of non-perturbative results beyond any 
finite number of loops in quantum field theory (QFT) 
perturbative expansion were obtained there. 

However, it should be emphasized that just because of 
the non-perturbative nature of the stochastic approach 
to inflation, it is based on a number of heuristic approxi- 
mations. Therefore, it is very important to check, when- 
ever possible, results obtained by its application using the 
standard perturbative QFT in curved space-time. Also, 
an inflationary space-time is not the exactly de Sitter 
one, H ^ 0, that can often lead to a drastic change in 



Description of the growth J2J in the free massless case in terms 
of the Fokker-Planck equation was first considered in Q. The 
stochastic approach was applied to the description of eternal in- 
flation in 
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conclusions. That is why our purpose here is to consider 
new applications of this approach to inflationary space- 
times with H 7^ 0, though with \H\ <§; H 2 (as required 
by observational data on a slope of the primordial power 
spectrum of scalar perturbations in the Universe) , and to 
compare results obtained in this way with those directly 
following from perturbative QFT in curved space-time. 
Our paper contains several novel results referring both to 
inflationary space-times and to the exact de Sitter space- 
time beyond one-loop approximation. Also, our results 
shed new light on the much discussed problem of the 
choice of an independent time variable in the Langevin 
and Fokker-Planck equations in the stochastic approach. 
Finally, we show that in most cases the mean square of 
a scalar field at last stages of inflation is very different 
from the instantaneous Bunch-Davies value ([T]) - even 
"eternal" inflation (which does occur in the inflationary 
models we consider for a sufficiently large initial value of 
an inflaton field) is not eternal enough for a light scalar 
field to reach equilibrium. 

The paper is organized as follows. In Sec. II the growth 
of fluctuations of a light test field during inflation driven 
by a massive inflaton is investigated. We prove that the 
results obtained in [l(| [ll[ using QFT with adiabatic reg- 
ularization in curved space-time can be obtained from a 
general diffusion equation with a noise term which has 
the same form as that in the Langevin equation for an 
inflaton field in the stochastic inflation approach (6(. In 
Sec. Ill more details are presented on how the mean 
square of a test scalar field with an arbitrary (though 
light, m -C H) mass is obtained using adiabatic regu- 
larization in curved space-time with a slowly changing 
curvature. Generation of light scalar field fluctuations 
in inflationary models can be so strong as to dominate 
finally the classical energy density of an inflaton and to 
drive a second stage of inflation. As an example, in Sec. 
IV we derive the conditions under which such a second 
stage can occur for the m 2 4> 2 inflationary model, with or 
without a break between the two stages of inflation. In 
Sec. V we discuss the impact of this inflationary particle 
production on the moduli problem. We show that quan- 
tum moduli problems are worse than the classical one 
and we improve on previous investigations (l2l [l3j . In 
Sec. VI we discuss the diffusion (Langevin) equation for 
inflaton fluctuations in case of a generic chaotic power- 
law potential V{4>) <x <fr n . Using the results obtained by 
QFT methods, we show that if one is interested in metric 
fluctuations, this diffusion equation should be formulated 
in terms of the independent time variable In a which is 
directly related to the number of e-folds during inflation 
N (this was among the possibilities envisioned in [(|). 
In Sec. VII it is shown that the stochastic method can 
also be used beyond the one-loop approximations for a 
self-interacting scalar field in the exact de Sitter space- 
time, again leading to the same results as obtained by 
QFT methods. We conclude in Sec. VIII and discuss the 
perturbative expansion of the diffusion equation in the 
Appendix. 



II. GROWTH OF LIGHT TEST FIELD 
FLUCTUATIONS DURING MASSIVE 
INFLATION IN STOCHASTIC APPROACH 

Let us first consider the production of particles and 
fluctuations for a test quantum field \ with a small mass 
m <C H. Thus, we neglect the x energy density and 
pressure in the background FRW equations. Also, it is 
assumed that there is no Bose-condensate of \, < x >— 
0. 

As an example, in this section we limit ourselves to the 
simplest case of inflaton with the quadratic potential. 
Some of QFT results for this model have already been 
obtained in [ljl [HI]. The Friedmann equation for the 
massive inflaton is: 
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where 8irG = M p 2 . In the slow-roll approximation, the 
inflationary trajectory is well approximated by 
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where the subscript denotes the beginning of inflation 
(these expressions were first presented in [14] in the con- 
text of a closed bouncing FRW universe). In the rest of 
the paper, we set ao = 1 for simplicity. The equation for 
inhomogeneous Fourier modes of \ is: 



Xk + 'AHxk 
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and we denote m 2 = a m 2 . Let us first discuss the case 
m x — m, or a = 1, in order to see the relation between 
renormalized QFT results and the stochastic approach. 
This case was studied in [l(| where inflaton fluctuations 
in a rigid space-time were considered. For the value of 
(x 2 ) renormalized by adiabatic subtraction, one obtains 
in the leading order: 
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This result agrees with the solution of the differential 
equation for (x 2 )ren: 
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(see e.g. Eq. (10) of paper [3j) with H not constant in 
time, but having the time evolution given by Eq. j4|). 
In the scope of the stochastic approach, the right-hand 
side of Eq. Q arises from a noise term in the Langevin 
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equation for the large-scale part of x, see Eq. ((36]) below, 
coarse-grained over the physical 3D-volume ~ (eH)~ 3 
which slowly expands during inflation. Here £ < 1 is 
an auxiliary parameter which enters into the definition 
of the coarse-graining scale fc cg = eaH in momentum 
space. Since we are interested in the leading order of the 
slow-roll approximation, terms of order H/H 2 coming 
either from taking an explicit time derivative of H or from 
solutions of the wave equation {7J at the course-graining 
scale may be neglected. With the same accuracy, one 
may alternatively take the physical course-graining scale 
k cg /a to be exactly constant during inflation and equal 
to eH(tf) where tf denotes the end of inflation. 

The general solution to Eq. ((9]) in the background (j4|) 
is indeed: 
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where C is an integration constant and we have used Eq. 
(6) in the second expression. For a long stage of inflation, 
when Hit) « H , we have [lOj]: 
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The regime H(t) << Hq occurs only when the number 
of e-folds is N w N = 3H 2 /{2m 2 ). 

The stochastic equation ([9]) can be easily extended to 
a generic m x not coinciding with the inflaton mass m: 
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Its solution is: 
(x 2 )ren = CiH 2 



C 2 H 2a 



3H{t) 



<x 2 > 



REN 



4tt 2 



(12) 



3H 4 



8ir 2 m 2 (a-2) 

3H 2a 
8TT 2 m 2 (2-a) 



{H^- 2a -H 4 - 2a )(13) 



where C\ , C% are integration constants. For the partic- 
ular marginal value a = 2, Eq. (I13|) should be replaced 
by: 
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When H(t) « H and a < 2, 
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In the limit a << 2 (but a ^ 0), \ renormalizcd fluc- 
tuations - and therefore the x energy density - depend 
non-trivially on the duration of inflation [10|, lllj . Note 
that only for a — 1 (i.e. m x — m) and at the beginning 



Hyi as occurs in the exact de Sitter 



of inflation, < x 
space-time. 

This shows the occurrenc e of a new characteristic scale 
in the problem, namely y/—H. In contrast to the de 
Sitter space-time, where H is the only scale present, in 



inflation we have also 



in the model un- 



der consideration. Fields which are light compared to 
the Hubble parameter, are subsequently divided in two 

classes: the ones with m x > y/— 3H and the ones with 

m x < y/—3H. Their corresponding particle production 
rates look different. 

A homogeneous solution of Eq. (T2")) given by the term 
containing C\ in Eq. (fTB")) is also the slow-roll solution 
for Xc\j with Xd being the classical homogeneous field 
x(t)- In the case a < 2, this part of the general solu- 
tion becomes negligible with respect to generated quan- 
tum fluctuations as inflation develops if (x 2 (^o))ren << 
3H$/ (87r 2 m 2 (2 - a)) (see Eq. CGl). In particular, if 
zero initially, (x 2 )ren becomes larger than the instanta- 
neous Bunch-Davies equilibrium value when 
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The energy density of x becomes comparable to that of 
the inflaton and cannot be neglected further on when 
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This never happens for a > 1, but for a < 1 back- 
reaction effects become important when 



H(t) < #(**) = H ( 
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Note that the dominant term in p x is rn 2 (x 2 )nEN, the 
X kinetic and gradient energy are sub-leading. Also, we 
always consider Hq < jM p \ with 7 = 0(1). 



It is important to note that (x 2 )ren is different 
from the instantaneous Bunch-Davies equilibrium value 
3iJ 4 /(87r 2 m 2 ) in Eq. for any a. Since H(t) decreases 
with time, it is only for a >> 2 that (x 2 )ren may ap- 
proach (but not become exactly equal to) the instanta- 
neous Bunch-Davies vacuum value at the last stage of 
inflation. 



III. RENORMALIZATION APPROACH 

We now wish to present the renormalized (x 2 )ren for a 
generic mass m x obtained using dimensional regulariza- 
tion and adiabatic subtraction. Such a calculation goes 
beyond [3, which addressed the first case in the previ- 
ous section (a = 1) only. 
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Let us consider solutions of the wave equation Eq. ([7]) 
which give time-dependent coefficients in the expansion 
of the quantum field \ m terms of the Fock annihilation 
and creation operators for each Fourier mode k, k = |k|. 
We use a set of approximate solutions of Eq. 0, chosen 
so as to separate the ultraviolet (UV henceforth) and 
infrared (IR henceforth) domains for k > eaH and k < 
eaH respectively (here e < 1 need not coincide with the 
parameter e in the definition of the coarse graining scale 
introduced above). 

The UV solution whose contribution will mostly be 
removed by the adiabatic regularization can be written 
in the generic form 



subtraction cancels the UV contribution more and more 
accurately as inflation develops. This defines a dynamical 
cutoff, and the IR part below this cut-off soon becomes 
dominant. Taking the time derivative of Eq. (f2"2"| . one 
obtains: 

^(0|x 2 |0W - ^ + 2a|<0| X 2 |0W • (23) 

This coincides with Eq. ([12]) for a generic mass to 2 = 
a to 2 since the direct computation of the integral in Eq. 
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for suitable functions v and A (whose values are normally 
close to 3/2 and 1 respectively, which are the values in 
the lowest order of the slow-roll approximation). 

The most important contribution is the far IR one. 
Such a solution can be computed in the same approxi- 
mation as used in our previous work [Tl| , where a match 
at the moment tk when k — ea(tk)H{tk) with the UV 
solution is imposed, together with the requirement of be- 
ing a solution of ((TJ in the deep IR limit k — ► 0. In the 
leading order, we choose the form 
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where the parameter x can be fixed using the second 
constraint. Note for the later use that at the moment tk 
one has 
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Substituting Eq. ([20]) into Eq. ||7J) and considering the 
deep IR region, one finds that x = 1 — a. 

In order to compute the renormalized value of (x 2 ) at 
the same scale as is used in the stochastic approach, we 
note that the contribution of IR modes should be taken 
up to k = ea(t)H(t) (the UV ones cancel due to adiabatic 
subtraction). Then it is natural to define tk by using e = e 
in Eq. (|2"Tj) in order to avoid terms like loge/e. We shall 
therefore make this choice. 

As mentioned above, the leading contribution to the 
renormalized value of (x 2 ) is given by the integration 
over IR modes. In the leading order of the slow-roll ap- 
proximation, one obtains the following result: 
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with I = eH$ and where the expansion of the Bessel func- 
tion for a small argument is performed. Thus, adiabatic 



in the slow-roll approximation leads to the second 
expression in Eq. ([T3"]) with C2 — 0. Therefore, the two 
methods of doing this calculation agree. 



IV. CONSEQUENCES OF PARTICLE 
PRODUCTION: SECOND STAGE OF INFLATION 

Let us discuss in more detail how the back-reaction 
of x may become important. If a > 1, then the back- 
reaction of x wm never be important during inflation 
(and afterwards, too, if the inflaton does not decay faster 
than x during (p)reheating). So, we consider a < 1 in 
the following. 

Let us introduce a new quantity @ = m 2 M^/H^. 
Back-reaction of x becomes important before or after the 
end of inflation driven by the inflaton <f> if H(t*) in Eq. 
([T8| is greater or smaller than to (which is the scale of 

) 2 inflationary model) re- 



tire end of inflation in the 



spectively. Therefore, x fluctuations become important 
before the end of inflation if 
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Thus, for a << 1 there is a second stage of inflation 
without a break if 
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and with a break if the < holds in the equation above. 
A second stage of inflation without a break starts in the 
presence of a light field with mass m x in the range 



1 >> a 
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which is a very broad range for large enough Nq. A 
second stage of inflation with a break would occur for 
smaller m x . In such a case the slow- roll of the inflaton 
field practically ends, but inflation would be sustained 
by quantum fluctuations of x corresponding to a much 
lower value for \H\ « m 2 /3 (see also [H,[iE[i3 f° r ex " 
amples of inflation as a whole, or its second stage driven 
by quantum fluctuations of a scalar field). Note that a 
second stage of inflation can occur even if curvature dur- 
ing the first stage of inflation is low, Hq << M p i, usually 
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after a break in this case. Thus, a second stage of in- 
flation seems to be very probable if the inflaton is not 
the lightest field around. Also, even without the second 
stage, x may play a role of a curvaton [H, [2(| Hl| . 



V. CONSEQUENCES OF PARTICLE 
PRODUCTION: IMPACT ON THE MODULI 
PROBLEM 

On considering x as a modulus, the above calculations 
allow a better estimate of the moduli abundance pro- 
duced by the expansion during inflation. The previous 
section shows that the moduli can even drive a second 
inflationary stage if a < 1 and the energy stored will 
dominate before the end of inflation (see Eq. (|18[) ). We 
focus here on the case for which \ does not dominate 
during inflation. 

In order to estimate the ratio n x /s (where n x = p x /m x 
and s are the number of \ particles and the entropy of 
produced particles, respectively), which is required to be 
less than 10~ 12 , we proceed as in We assume the 

immediate thermalization for the inflaton and x after the 
accelerated expansion: 



2n 2 gT 3 
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where T is the reheating temperature and g is the number 
of species (for instance g = 106.75 for the standard model 
d). 

We first consider a >> 2 and obtain: 
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Because of strong positive dependence on the reheating 
temperature, the requirement n x /s < 10 -12 can easily 
be satisfied. 

We now consider a < 2, the case for which particle 
production is very different from the one extrapolated 
from the exact de Sitter space-time. The energy density 
stored in x f° r H << Hq is: 
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This equation replaces Eq. (4.16) of [ll 



(30) 

t Note that 



t Eq. (4.16) of |l.'J agrees with Eq. (130 I t for a = 1 only, neverthe- 
less, it is used for any a. Also, the growth law JJJ is used there 
which is valid for a = 1 only as pointed above. 



the dependence on T is rather peculiar: for a > 3/4, it 
grows with T, while for a < 3/4, it decreases with T. For 
a ~ 0, the ratio is 
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Therefore the problem of light moduli cannot be solved 
by lowering the reheating temperature in the m 2 4> 2 model 
and is much worse than expected in [l3l |. 



VI. INFLATON FLUCTUATIONS AND 
STOCHASTIC APPROACH FOR GAUGE 
INVARIANT FLUCTUATIONS 



An inflaton has an effective mass which is much smaller 
than the Hubble parameter during inflation. Therefore, 
the above approach should hold for inflaton fluctuations, 
too. However, such a case differs from the previous one, 
since scalar metric fluctuations should be taken into ac- 
count in addition to field (inflaton) ones. In other words, 
the quantity which is quantized is a linear combination 
of field and metric fluctuations. As in [fl|, we choose 
the gauge in which inflaton fluctuations coincide with 
the gauge-invariant Mukhanov variable which is canoni- 
cally quantized 23] . Then our results will be valid in any 
gauge, if a gauge invariant variable is considered. 

The time evolution of the renormalized inflaton fluctu- 
ations for the m 2 (j) 2 inflationary model has been already 
obtained in [III, [2J] using a perturbative QFT analysis of 
the Einstein equations. To the lowest order in the slow- 
roll and in long- wavelen gth approximations, equations for 
fluctuations in the first [ll[ and second [2J| order are: 



6H 



(32) 



6H 



2HM 2 l 



(V 1} ) 2 ■ (33) 



We now observe that these two equations can be ob- 
tained order by order from the expansion around a clas- 
sical solution, (f>(t, x) = 4> c i (t) + S^ 1 ' + Sc/)^ , of the equa- 
tion: 



dN 
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where N = log(a/ao) is the e-fold number from the be- 
ginning of inflation. This suggests that the equation (|34|) 
is valid to all orders in S(f> = <p — <t>ci (t) in the context of 
the approximation used. Returning to the time evolution 
and using dN = Hdt, one has: 
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where = and so on. The quadratic inflaton po- 
tential leads to the result in Eqs. ((32)) and (|33|) after 
expanding 5<j) = 5<jP-' + S(f>^ and taking the coefficients 
in the brackets to the lowest order in slow-roll. The same 
method may be used for an arbitrary inflaton potential 
V(cf>). 

How does one re-derive these results in the stochastic 
approach? The consideration above suggests that one has 
to choose the time variable N = J H(t)dt in the Langevin 
stochastic equation for the large-scale part of <f> for this 
aim. Then it acquires the form (see Eq. (55) of @): 
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where H 2 (t) = V(0(f))/3M*, in this (leading) order of 
the slow-roll approximation, thus, H may be considered 
as a function of </>. As proved in [6| (see also [9[), though 
formally the noise f(N) is an operator quantity contain- 
ing the Fock annihilation and creation operators ak and 
at for modes with k = eaH, its values in different points 
of space and for different N are commutative. Thus, 
it is equivalent to some classical noise whose distribution 
function appears to be Gaussian since the quantum noise 
/ is linear in ak and a\ . Correspondingly, the large-scale 
part of the quantum inflaton field <f> is equivalent to a 
classical stochastic field <fi with some normalized proba- 
bility distribution p((f>, N), J p((f>, N) d<p = 1, in the sense 
of equality of all possible classical and quantum expecta- 
tion values: 



REN 



(37) 



where F means any functional which may include time 
and spatial derivatives of any order. * 

Following 6] , the number of e-folds N was considered 
as a time variable in the stochastic Langevin equation 



t Note that the so called volume- weighted averaging |25j , for which 
p is not normalizable and thus loses the sense of a probability 



in a number of papers, e.g. in (28| and most recently 
in [2{|, while in many other ones the proper time t was 
used, e.g. in [25|, [30, bLU, [32| . This usage of different time 
variables should not be mixed with the invariance of all 
physical results with respect to a (deterministic) time 
reparametrization t — > f(t) which is trivially satisfied af- 
ter taking the corresponding change in the metric lapse 
function into account. In contrast, the transformation 
from t to N made using the stochastic function H(<j>(t)) 
leads to a physically different stochastic process with an- 
other probability distribution. Our new statement in this 
paper following from the arguments above is that one 
should use the N variable when calculating mean squares 
of any quantity containing metric fluctuations like the 
Mukhanov variable 5(f) or the gauge-invariant metric per- 
turbation £. Otherwise, incorrect results would be ob- 
tained using the stochastic approach which would then 
not coincide with those obtained using perturbative QFT 
methods. This statement is further supported by exact 
non-perturbative results (valid to all orders of metric per- 
turbations) from the general dN formalism which relates 
the value of £ after inflation to the difference in the num- 
ber of e-folds N in different points of space - it was first 
used in Q for one inflaton field, see Eq. (17) of this pa- 
per, and then generalized to multicomponent inflation in 
[H, HI], see also [35|,[3(| for its recent developments. 

On the other hand, the usage of the t variable is nat- 
ural if one is interested e.g. in differences of the local 
duration of inflation in proper time in different points 
of space which, in principle, may be measured using a 
spatial distribution of a phase of the wave function of a 
heavy particle with a mass m H. Thus, the choice of 
a proper time variable in the stochastic equation is not 
an absolute one but is dictated by the physical nature of 
'clocks' relevant to observable effects. In our case, N is 
the 'clock'. E.g., if a rather small (~ 10%) contribution 
from the integrated Sachs- Wolfe effect due to the exis- 
tence of dark energy (or a cosmological constant) in the 
present Universe is neglected, a large angle anisotropy of 
the cosmic microwave background temperature for mul- 
tipoles / < 50 is just given by AN = N — (N) at the last 
scattering surface: AT/T = —AN/5, so that a larger lo- 
cal amount of inflation produces a negative spot in the 
CMB temperature map (here () means averaging over the 
2-sphere - the last scattering surface). 

Now, in order to compare with perturbative QFT 
results, let us assume that fluctuations in <f> are still 
much less than the classical background value <fi c i(N): 
\6<fi\ = \4> — <f>ci\ "C \4>d\- In the rest of this section, we 
shall denote by H and <j> their classical background val- 
ues. Then the noise term in Eq. (|36p may be considered 



distribution, is nowhere used in this paper, as well as in 
For this reason, criticism of this approach contained e.g. in the 
recent paper [26| has no application to our results. See [27] for the 
most elaborate proposal of how to cure problems of the volume- 
weighted approach. 
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as a perturbation. In the first order, after expanding the 
first term in the right-hand side of Eq. (|36p in powers 
of 5(f), we get the following equation for 5 = 5<jft~> (the 
superscript 1 denotes the order of expansion) : 



(38) 



a family of chaotic inflationary models with potentials 
V{4>) — —K n ~ 4 (f) n and obtained for such models results 
similar to those presented in section II, but valid for 
gauge invariant fluctuations. Eq. (|40[) has been found 
to give renormalized values in complete agrement with 
our perturbative QFT calculations in the slow roll ap- 
proximation. 



where the prime denotes the derivative with respect to 
4>. Multiplying both sides of Eq. ([38]) by 6, averaging 
and using the relation (f5) = H 3 /8ir 2 , we obtain the 
equation for u = {(S(j)^) 2 } corresponding to the one-loop 
approximation of QFT in curved space-time: 



du 
dN 



H 



H 2 

= (39) 



Eq. (f3"5)) is valid for any potential V{<f>) satisfying the 
slow- roll conditions. Its generic solution is 



H 



12 



%-n 2 B 2 Ml l 



(40) 



where the integration variable has been changed from N 
to <j) using the slow-roll relation dN = -H/(2H'M 2 l ) d<\>. 
Applying this solution to the case V(<p) = m 2 (j) 2 /2, H — 
m<f> I ' (sf&Mpi) and assuming u(0) = 0, we obtain: 



it = 



M 6 



H 6 



8ir 2 m 2 H 2 



(41) 



that just coincides with the result derived in [llj using 
QFT methods. Also, applying Eq. (|40)l to power-law 
inflation where a(t) ~ t p with p » 1 and H exponen- 
tially depends on 0, we easily recover the correct leading 
renormalized value of u obtained in [37l |. 

In the Appendix, the mean value of a second order 
field fluctuation (Scj)^) is calculated and also conditions 
of the validity of the perturbation theory are reviewed. 

To emphasize the difference, note that would we re- 
peat the same procedure using the stochastic equation 
(I36p written in terms of the independent variable t, 
dt = dN/H, we would obtain a different result |25l. 1321]: 



H r 



8ir 2 M 2 pl 



H ,3 



which, in particular, reduces to 3(Hq — H 4 ) / (16ir 2 m 2 ) 
for massive inflation. As already explained above, the 
discrepancy between (|40]) and (f42|) is not surprising and 
reflects the fact that different stochastic processes are 
considered, that is why we wrote u instead of u in Eq. 

The general result presented in Eq. (|4"0"]) for the renor- 
malized value of the average squared gauge invariant 
massive inflaton fluctuation has been explicitely veri- 
fied by the authors in the context of QFT with renor- 
malization obtained on employing the adiabatic sub- 
tration prescription. In particular we have considered 



VII. 4-LOOP CALCULATION FOR A 
MASSLESS SELF-INTERACTING TEST SCALAR 
FIELD IN THE STOCHASTIC APPROACH 

In the previous sections, only the mean value of 
{(54>^) 2 ) corresponding to one scalar loop in external 
background curved space-time was calculated (with small 
metric fluctuations taken into account also for the case of 
an inflaton scalar field). However, as was shown in 0, 
the stochastic approach can reproduce QFT results for 
any finite number of scalar loops and even beyond (e.g. 
results obtained using instantons). As an example, let 
us consider a massless self-interacting test scalar field \ 
with the potential V(\) = A% 4 /4 in the exact de Sitter 
space-time with the curvature Hq following [9(. Then 
N = Hot, and it makes no difference which time vari- 
able is used in the stochastic Langevin equation (|36|) . It 
is straightforward to construct the corresponding Fokker- 
Planck equation for a normalized probability distribution 



dp 
9iV 



A d( PX 3 ) 



H 2 d 2 p 



(43) 



3F 2 d x Sir 2 d x 2 ' 

Multiplying Eq. (|43[) by where n is an even integer 
and integrating over x from — oo to oo, we obtain the 
following recurrence relation (c.f. also [381]): 



dN 



(x n ) 



3Hi 



(x n+2 ) + 



n(n - 1)H 2 
8^2 



(X" 



(44) 



Then solving Eq. (I44[) iteratively beginning from n — 2 
with the initial conditions (x"(0)) = for all n, we find: 



(x 2 



'REN 



Air 2 { 



aiXX 2 + a 2 \ 2 X i + a 3 A 3 A e 



(42) <X 4 )ren = 3 (^pj (l + p 1 \X 2 +(3 2 \ 2 X i + ..) 



(x 6 



'REN 



= 15 



Air 2 



(I + 71AA 2 + ...) 



(45) 



where X = N/(2ir) and 
2 

ai = "3 
Pi = -2 0a 



a 2 



4 

" 5 
212 

l5~ 



a 3 



71 



424 
315 

-4. 



(46) 



This one point distribution function should not be confused with 
the energy density p x of a test \ field considered in Sects. II-V. 
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The result for ai a gree s with the perturbative QFT 

computations in [33, (see also [4l[). Since the A- 
independent term in the expression for (x 2 )ren corre- 
sponds to one scalar loop and each next power of A re- 
quires one more scalar loop, the total result for (x 2 )ren 
presented in (145|) and (|46|) requires consideration of di- 
agrams up to 4 loops included in the perturbative QFT 
treatment. 



VIII. CONCLUSIONS 

We have studied the application of the stochastic ap- 
proach to inflationary space-times with H ^ and 
compared its results to those obtained using the stan- 
dard QFT in curved space-time, extending the results 



of [H lllj We can summarize our main results by the 
following points: 

1. On considering a test field \ with mass m x in 
m 2 (jr inflation, we have shown that the instantaneous 
Bunch-Davies expectation value (\ 2 ) = 3ii 4 (f)/(87r 2 m 2 ) 
is never reached. It may be approached at the end of 
inflation for m x >> m only. If m x << m, the value of 
(x 2 ) at the end of inflation is quite different from that 
extrapolated from the exact de Sitter space-time. 

2. We have analysed implications of the particle pro- 
duction peculiar to m 2 (j) 2 inflation. The moduli problem 
is more serious than in the classical counterpart and also 
with respect to previous quantum investigations [l3|. 

3. Concerning gauge-invariant inflaton fluctuations, 
we have clarified why the stochastic Langevin equation 
for the large-scale part of an inflaton field should be for- 
mulated using the number of e-folds N — In (a/do) as 



an independent time variable, if one is interested in any 
result regarding gauge-invariant inflaton fluctuations and 
metric perturbations. The mean square of inflaton fluctu- 
ations calculated in this way has been shown to coincide 
with the earlier result of obtained using standard 
perturbative methods. 

4. In the inflationary space-times studied here, neither 
test fields nor inflaton mean square expectation value ad- 
mit a static equilibrium solution. 

5. The equivalence between the stochastic and the 
standard field-theoretic approaches works beyond the 
one-loop approximation to QFT in curved space-time. 



Acknowledgements 

AS was partially supported by the Russian Foundation 
for Basic Research, grant 08-02-00923, and by the Re- 
search Programme "Elementary Particles" of the Russian 
Academy of Sciences. This work was started during visits 
of AS to Bologna in 2003-5 financed by INFN: we thank 
INFN for support. Some of these results were presented 
at the International Conference on Theoretical Physics in 
the Physical Lebedev Institute, Moscow (Russia), April 
11-16, 2005. 

IX. APPENDIX 

In this appendix, we shall denote by H and <f> their 
classical background values; we also use u = {{8(f>^) 2 ). 
Starting from Eq. (|36|) , and performing a perturbative ex- 
pansion, it is easy to derive for the average second order 
fluctuation (S(j)^ ) a time evolution also governed by (see 
also [H) 



d_ 



2)\ 



H 3 
16tt 2 



3H 



4"t> 



H 
''H 



2)\ 



3H 



-,V A . 



H 



-.Vs. 



H 



V 



(47) 



The general solution of Eq. (|4T)) with the initial condition 8<f>(t =ti) = is given by 



dt' 



16tt 2 



3H 



-.Vs. 



H 



-.Vs. 



H ) V 



(48) 



In the particular case of a chaotic inflation with V(<f>) — 
2^-0 2 we obtain Eq. (T4"Tj) for first order and: 



Now we wish study the validity of the perturbative 
expansion by considering the ratio ^ and . Using 
the former results one obtains 



(2)\ _ 



32Tr 2 m 2 HM 2 l 



Hi 



H 2 



-3K-H 4 ) 



(49) 

for second order. Note that the solution (T4"T]) clearly 
agrees with the solution (fT3"|) for suitable initial condi- 
tions and ot = —1. 



2)\ 



1 



1 ^(H$-H e 



8ttV3 M, 



1 



P i 



1 1 



(#o 6 



if 6 ) 



1/2 



4ttV3M p1 m 



{Hi - H 



(50) 
(51) 
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Thus, as we can see, those two ratios differ only by a 
factor 2 for the leading term toward the end of inflation. 
To see when the perturbative expansion is no longer use- 
ful we can use the variable N defined as the number of 
c-folds away from the maximum value N max = No = 

!°gfl77y = iS ^mely 



N = log 



oft) 
a(U) 



2 m 2 



(52) 



If we consider the end of inflation when H = m this cor- 
responds to N = 3/2, so in order to have a perturbative 
expansion with small terms for all of the duration of in- 
flation we obtain, on considering the condition (|54")) with 
N = 3/2, the following condition on Hq: 



/24tt 



H < [ -^Mptm 2 



1/3 



(55) 



Using this variable one obtains, to leading order, for the 
ratio (J3UJ) the following result 



2) 



) y/2 m N„ 



24tt M„i N 



(N max - N) 



~\l/2 

AM 



(53) 



If we require that the absolute value of this ratio be less 
then one we obtain, under the condition 7687t 2 -t# >> 1 

"0 

the following approximate constraint 



TV > 



3V3 



n o 



48 7T Mpim 2 



(54) 



similarly on considering the ratio (j5"Tj) one obtains the 
condition 



H < \-y=M p im' 



1/3 



(56) 



the above correspond, for the particular numerical value 
M p i = 10 5 m, respectively to Ho < 163.28m and Hq < 
129.596m. Such bounds are in agreement with previous 
investigations 24] . 
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